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ABSTRACT 

The census of young moving groups in the solar neighborhood is significantly in- 
complete in the low-mass regime. We have developed a new selection process to 
find these missing members based on the GALEX All-Sky Imaging Survey (AIS). 
For stars with spectral types >K5 (R — J > 1.5) and younger than ?^300 Myr, 
we show that near-UV (NUV) and far-UV (FUV) emission is greatly enhanced 
above the quiescent photosphere, analogous to the enhanced X-ray emission of 
young low-mass stars seen by ROSAT but detectable to much larger distances 
with GALEX. By combining GALEX data with optical (HST Guide Star Cat- 
alog) and near-IR (2MASS) photometry, we identified an initial sample of 34 
young M dwarf candidates in a 1000 sq. deg. region around the ^10-Myr TW 
Hydra Association (TWA). Low- resolution spectroscopy of 30 of these found 16 
which had Ha in emission, which were then followed-up at high resolution to 
search for spectroscopic evidence of youth and to measure their radial velocities. 
Four objects have low surface gravities, photometric distances and space motions 
consistent with TWA, but the non-detection of Li indicates they may be too old 
to belong to this moving group. One object (M3.5, 93±19 pc) appears to be the 
first known accreting low-mass member of the f=al5 Myr Lower Centaurus Crux 
OB association. Two objects exhibit all the characteristics of the known TWA 
members, and thus we designate them as TWA 31 (M4.2, 110±11 pc) and TWA 
32 (M6.3, 53±5 pc). TWA 31 shows extremely broad (447 km s _1 ) Ha emission, 
making it the sixth member of TWA found to have ongoing accretion. TWA 32 is 
resolved into a 0.6" binary in Keck laser guide star adaptive optics imaging. Our 
search should be sensitive down to spectral types of at least M4-M5 in TWA and 
thus the small numbers of new member is puzzling. This might indicate TWA 
has an atypical mass function or that the presence of lithium absorption may be 
too restrictive a criteria for selecting young low-mass stars. 

Subject headings: Stars: activity, chromospheres, coronae, late-type, ages - Sur- 
veys: Ultraviolet - Galaxy: solar neighborhood 

1. Introduction 



Observational studies of planet formation ha ye been energ i zed by the discove ry of young 
(<100 Myr) solar-type stars close to Earth (e.g. I Jeffries! Il995t IWebb et al.lll999l ). identified 
from multiple indicators of youth including chromospheric activity and strong X-ray emission. 
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The combination of distances, proper motions, and radial velocities (RV s) has allowed many 



of th e se stars to be kine matically linked to coeval moving groups (e.g., iZuckerman &: Song 
2004 ; iTorres et al.l 120081 ) . These young moving groups (YMGs) are several times closer to 
Earth than the traditional well-studied star-forming regions such as Taurus and Orion (~150- 
500 pc). More importantly, these groups have ages of ~10-100 Myr, a time period in stellar 
evolution that has been largely underrepresented in previous studies. This is expected to 
be a key epoch for understanding planet formation, coinciding with th e end of giant planet 
formation and the active phase of terrestrial planet formation (e.g. iMandell et al.l 120071 ; 
Ida fc Linllioosh . 



One of the best studied of these YMGs is the TW Hydrae Association (TWA;lKastner et al 



1997). Its unique combination of youth (~10 Myr; iBarrado Y Navascueal2006l ; lMentuch et al 



20081 ). proximity (~30-100 pc) and (relative) compactness on the sky, makes it is a partic- 
ularly promising observational test bed for theories describing disk evolution and planet 
formation. T WA is named for the actively accreting K7 star with a circumstellar disk 



(jHenizd 1 19761 ; iRucinski fc Krautterl Il983l ; Ide la Reza et al.l Il989l ) . The discovery of such a 
young star far from any star-forming region prompted several searches for additional T 
Tauri stars in the same par t of the sky. These searc hes used infrared excesses with IRA S 
jGregorio-Hetem et al.lll992h. ROSAT X-rav activity Jsterzik et al.l ll999l: IWebb et al.lll999h , 



and common space motion dSong et al. 



2003 



Scholz et a 



20051 ). plus near-IR photometric 



searches for brown dwarfs (IGizisI 120021 ; lLooper et al.l 120071 ) . The current census of accepted 
TWA members is up to 25 objects (plus 2 discovered by us), including either late-K or 
early-M stars and a few brown dwarfs, making TWA seem unusua l in its initial mass 



tion compared to other star forming regions (e.g. ITorres et al.ll2008l ; ISlesnick et al.ll2008) . In 



particular, there appears to be a dearth of M3-M7 stars (Figures [T] and [21 iMamajek 



unc- 



2005 



and references therein). It is possible that this large incompleteness arises from the fact that 
most young star searches have mostly relied on bright optical catalogs (e.g., Hipparcos) as a 
starting point, which are deficient in M dwarfs, or perhaps TWA may have an unexpected 
mass function. 

Over the last few years, we have been undertaking a systematic effort to complete the 
young low-mass census through multi-catalog selection and followup spectroscopy. Our initial 
work focused on the immediate solar neighborhood (<25 pc), identifying the youngest M 
dwarfs from the known (Gliese & NStars) catalogs. Our spectroscopy had a very productive 
Vo) confirmation rate, identifying 144 new young (< 30 Myr) M dwa rfs using X-ray 



select ion (ROSAT Brig ht Source Catalog; IVoges et al.lll999l ) as a first cut (jShkolnik et al. 



20091 . hereafter SLR09; IShkolnik et al.ll2008l . BOlOh . By scrutinizing the GALEX NUV and 
FUV properties of this 25-pc sample, we have constructed efficient selection criteria to expand 
the search for young M dwarfs to 100 pc thereby making the M dwarfs of many of the YMGs 
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accessible. In this paper, we apply our new methodology to TWA to find its 'missing' M 
stars. 



2. Selection of Young Stars using the GALEX All-sky Imaging Survey 



Stellar activity is a well-established indicator of youth (e.g.. |Preibisch fe Feigelsonll2005l ). 
and the only all-sky surveys suited for finding active, young stars have been the ROSAT X- 
ray catalogs (e.g. Voges et al. 1999). However, since the luminosities of M dwarfs are ~10- 
300 x lower than solar-type stars, ROSAT is generally limited to the nearest, earliest-type 
M dwarfs. 



The NASA Galaxy Evolution Explorer (GALEX; iMartin et al.l 120051 ) provides a new 
resource that enables a major expansion of the young low-mass census, far beyond previous 
data sets. The GALEX satellite was launched on April 28, 2003 and has imaged most of the 
sky simultaneously in two bands: near-UV (NUV) 1750-2750 A and far-UV (FUV) 1350- 
1750 A, with angular resolutions of 5"and 6.5" acr oss a 1.25° field o f view. The full description 
of the instrumental performance is presented by iMorrissey et al.l (120051 ). The GALEX mis- 
sion is producing an All-sky Imaging Survey (AIS) which is archived at the Multi-mission 
Archive at the Space Telescope Science Institute (MAST)Q The NUV/FUV fluxes and mag- 
nitudes are produced by the standard G ALEX Data Analysis P ipeline (ver. 4.0) operated at 
the Caltech Science Operations Center ( Morrissey et al.l I2OO5I ) . The data presented in this 
paper made use of the fourth data release of the AIS (GR4), covering 2/3 of the sky@ 

For M dwarfs, the flux in the GALEX bandpasses is m ade up of an abundan ce of 
emission lines. St ellar flare activity on M dwarfs observed by iRobinson et al.l (120051 ) and 
Welsh et al.l ( 120071 ) has shown that GALEXs FUV flux is mainly due to transition region 
C IV (~50%), while weaker emission lines and continuum provide the remaining 50%. In 
the NUV band, flux is primarily due to continuum emission plus Mg II (10%), Fe II (17%), 
Al II and C III (14%) lines. I n addition, co ronal lines provide 2% of the FUV and 10% 
of the NUV emission. (See alsoEaganoll2009h Th is makes GALEX ideal for finding active 
low- mass stars (e.g. iFindeisen fe Hillenbrand! l2010h . 



We first tested GALEXs sensitivity for finding young M dwarfs by correlating the X- 
ray-selected sample of M dwarfs within 25 pc described in SLR09 with the GALEX/ AIS 



1 One can query the AIS through either CasJobs ( http: / /mastweb. stsci.edu/gcasjobs J\ or a web tool called 
GalexView (http://galex.stsci.edu/galexview/). 

2 See details at http://www.galex.caltech.edu/researcher/techdoc-ch2.html 
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archive. Sixty-seven percent were detected in the NUV band, corresponding to the 67% 
sky coverage of the AIS in the GR4 data release, i.e. all those observed by GALEX were 
detected. This indicates that GALEX is at least as sensitive as ROSAT in identifying young 
M dwarfs. The GR4 release also detected 17 of the 25 (68%) known M dwarfs in TWA (not 
double-counting visual binaries (VBs) not resolved by GALEX), with 7 of the 8 remaining 
not yet observed. The o nly TWA members ob served by GALEX but not confidently detected 
were TWA 30 A and B Looper et al.lhoiOal lM 



We then carried out a search for UV detections with the NStars 25-pc census (~1500 M 
dwarfs; Reid et al. 2007). Those with NUV detections are plotted in Figure [3] as a function 
of R — J color as a proxy for effective temperature. From the results of this GALEX 25-pc 
analysis, we found that NUV data yield many candidates, but applying FUV criteria provide 
an excellent means to distinguish between the (never-before-delineated) quiescent emission 
of old stars (Fpuv / ' Fj < 10~ 5 ), the faint sources (FUV not detected), and the truly young 
targets with high levels of NUV and FUV emission, where F^ruv, Fpuv and Fj are the 
stellar fluxes in the respective bandpasses. We find that for R — J <4 (SpT < M4), we can 
detect young M dwarfs at least out to several hundred parsecs. And for young stars later 
than R — J«6 (SpT ~ M7), we can probe out to distances of ~75 - 100 pc and further for 
very active stars. Thus GALEX offers an enormous advantage for detecting young M dwarfs 
compared to ROSAT. 

In principle, flaring M dwarfs could contaminate our sample since UV flares would be 
mistaken for strong steady-state UV excesses. However, we expect the contribution of these 
objects to be less than 3%, based on the fraction of short-term flaring M dwarfs from GALEX 
results thus far (Welsh et al. 2006). 

To summarize, the GALEX/ AIS is substantially more sensitive than previous ROSAT 
catalogs, meaning that young M dwarfs can be identified at much greater distances and to 
later spectral types. In this study, we use the GALEX NUV and FUV data combined with 
existing optical and NIR photometric catalogs to identify TWA candidate members within 
100 pc of the Sun@ 



3 GALEX does have a weak (2.4er) detection in the NUV bandpass 26" away from TWA 30 A's 2MASS 
coordinates, with no corresponding FUV detection. The non-detection of TWA 30 A may be related to the 
unusually low Ha EW and spectral variability observed for this target. lLooper et al.1 (|2010aj) speculate this 
is due to an accretion disk viewed nearly edge-on with the stellar rotation a xis inclined to the dis k. Similarly, 
there is an 8a detection 27" away in a different direction from TWA 30 B (jLooper et al.ll2010bl ). 



4 While this manuscript was under review, a preprint bv lRodriguez et al.l (|2010h presented a NUV+NIR 
search for members of TWA and the Scorpius-Centaurus, including the star we refer to as TWA 32. However, 
their RV and UVW of this star are inconsistent with our measurements. They exclude it as a member of 
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3. Sample Selection of TWA Candidates 



Late-K and M dwarfs have distinctive photometric properties, and we can use these to 
identify stars within 100 pc. We combin ed optical photometry from the HST Guide Star 
Catalogue (GSC 2.3: lLasker et all l2008h with JHK photometry from the 2MASS Point 
Source Catalogue (ICutri et al.ll2003l ) to identify candidate late-type dwarfs. We first queried 
the 2MASS catalog for all sources with H — K > 0.25 (FigureEJ, aiming to include everything 
with SpT later than M2 within a 1000 sq. deg. region around the known TWA members 
bounded by these position limits: RA=10 to 13 hrs, DEC=-30 to -60 degrees, b > 10 
degrees (Figure 0J. This yielded 261,547 objects. Cross-matching these against the GSC 
2.3 catalog returned 183,361 targets. A 10" cross-matching with the GALEX/AIS returned 
1968 targets with >3a detections in both the NUV and FUV, as well as eliminated any 
background early-M giants. 

The distance li mit was set using an optical/near- infrared color-magnitude diagram, Rp 
vs. Rp-K (Figure 151 iReid et al.ll2007l ). where Rp is GSC 2.3's photometric optical band with 
a wavelength range from 6000 to 7500 A, very similar to the standard R passband. This 
method has proven extremely effective at identifying cool dwarfs within 25 pc, and it was 
a simple matter to extend it to larger distances. For young pre-main-sequence (PMS) stars 
which are over-luminous compared to dwarfs, the effective distance limit becomes 120 pc for 
the early M's based on PMS models by lBaraffe et al.l ( 119981 ) . Our resulting photometric cuts 
were: 



R F < 2.85( J R F -A)+5.94; 2 < (Rp-K) < 3.7 
R F < 6.5(i? F -A)-7.55; 3.7 < (Rp-K) < 3.9 
Rp < 1.76(^-/0+10.94; 3.9 < (Rp-K) < 6.5 

Fifty-five targets remained after applying these color-cuts. The TWA candidates chosen 
have both high NUV and FUV fractional luminosities (Fpuv / Fj > 10~ 5 and F^uv/Fj > 
10~ 4 ), comparable to or greater than strong X-ray emitting M dwarfs (SLR09) and known 
TWA members. This generated a target list of 38 stars (including 4 previously known TWA 
members, TWA 3A, 3B, 10 and 12)@ for spectroscopic followup (Figures [3] and |5J. 



any known YMG. 

5 The known TWA members not listed here were outside of our search criteria, i.e. had H — K < 0.25, 
had declinations above -30 degrees, or, as is the case for the brown dwarf members, were not detected in 



-7- 



4. Spectroscopic Confirmation 



We then subjected our compiled list of GALEX-selected TWA candidates through a two- 
step process of ground-based spectroscopy to assess their potential membership in TWA. We 
first acquired low-resolution spectra to search for th e activity-sensitiv e emission line Ha 
always present in young low-mass stars (<400 Myr; IWest et al.l 120081 ). T hirty of the 34 



cand idates were observed with the Magellan Echellette Spectrograph (MagE; Marshall et al. 



20081 ) mounted on the 6.5-m Clay telescope on UT 2009 February 10 (Table [T]). MagE is a 
moderate resolution, cross-dispersed echellette, covering the optical wavelengths from 3100 - 
10000 A. We used the 0."7 slit which produced a resolution of 4100. Wavelength calibration 
was performed with the ThAr exposures taken before and after each stellar target. 

Of these 30, 16 exhibited Ha in emission (EW <— 1 A), for which we acquired high- 
resolution spectra using the Magellan Inamori Kyocera Echelle (MIKE) spectrograph also at 
the Clay telescope over several nights: UT 2009 April 15, June 6-8, December 31 and 2010 
January 1. We also observed 14 of the known TWA M dwarfs^] (Tables [2] and |3]). We used the 
0.5" slit which produces a spectral resolution of ~35000 across t he 4900 - 10 000 A range of 



the red chip. These data were reduced using the facility pipeline (jKelsonl 120031 ). Each stellar 
exposure was bias-subtracted and flat-fielded for pixel-to-pixel sensitivity variations. After 
optimal extraction, the 1-D spectra were wavelength calibrated with a ThAr arc. To correct 
for instrumental drifts, we used the telluric molecular oxygen A band (from 7620 - 7660 A) 
which aligns the MIKE spectra to 40 m s _1 , after which we corrected for the heliocentric 
motion of the Earth. The final spectra are of moderate S/N 25 per pixel at 8000 A). 

The MIKE data provide RV measurements to better than 1 km s _1 in almost all cases. 
We measured RVs of our 16 candidates and 14 known TWA members by cross-correlating 
their spectr a with those of tw o RV standards taken on t he same night, nam ely GJ 179 
(SpT = M4; iMarcv et al.lll987h and/or GJ 908 (SpT = Ml; iNidever et al"ll2002h . We cross- 
correlated each of 9 spectral orders between 7000 and 9000 A (excluding those orders with 
strong telluric absorption) whe re M dwarfs emit most of their optical light using IRAF'fl 
fxcor routine (lFitzpatricklll993l ). We measured the RVs from the gaussian peak fitted to the 
cross-correlation function of each order and adopted the average RV of all orders. The target 
RVs and their standard deviations are listed in Table HJ 



the FUV bandpass at the 3ct level. 

6 Our spectrum of TWA 23 revealed it to be a double-lined spectroscopic binary (SB2). 

7 IRAF (Image Reduction and Analysis Facility) is distributed by the National Optical Astronomy Ob- 
servatories, which is operated by the Association of Universities for Research in Astronomy, Inc. (AURA) 
under cooperative agreement with the National Science Foundation. 
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The high-resolution spectra also allowed us to identify multi-lined binaries whose excess 
UV emission may be due to the tidal locking of a close-in binary system rather than youth. 
We found 2 (possibly 3) of the 16 to be SB2s con sistent with the 16% l ow-mass spectroscopic 
binary fraction of our ROSAT study (SLR09, IShkolnik et al.l |2010| ). With single epoch 
observations we are unable to exclude SBls from the remaining sample, however, given the 
results of SLR09, we do not expect any to be found with orbital periods less than 1-2 years. 



5. Spectral Types 



The spectra of M dwarfs are dominated by the strong TiO molecular bands, which 
are particularly diagnostic of the sta r's temperatur e . To estimate spectral types (SpT) we 
used the TiO-7140 index defined by IWilking et al.l (120051 ) as the ratio of the mean flux in 
two 50- A bands: the 'continuum' band centered on 7035 A and the TiO band on 7140 A. 
We calibrated the TiO-7140 in dex for our data sets using 126 M dwarfs from the NStars 
25-pc sample (IReid et al.l 120071 ) we observed with MIKE that have published spectral types 



plus several RV sta ndards and known members of TWA ( IReid et al.l Il995l ; iRiaz et al.l 12006 



Wenger et al.ll2007l and references therein). The relationship we used to convert the TiO-7140 
index to SpT for M0 - M6 stars is: 



SpT=(TiO 7 i4 - 1.1158)/0.1802, rms= 0.5 subclasses (1) 



Here, M0 corresponds to SpT=0, Ml — > 1, M2 — > 2, etc. The linear fit to the calibra- 
tion is consistent to better than 0.2 subclasses with that used for the Keck/HIRES and 
CFHT/ESPaDoNS data presented in SLR09. We determined the errors of our measurements 
for the TiO index (and subsequent indices discussed below) by taking the mean standard de- 
viations measured for 5 RV standards stars observed with the same setup on multiple nights. 
We calculated an average error of 0.012 in the TiO index measurement, which translates to 
only 0.07 subclasses in SpT. Although this uncertainty is small, the calibration is based on 
a sample with SpTs binned to half a subclass, imposing a 0.5 subclass uncertainty in the 
calculated SpTs listed in Table [TJ 



6. Identifying the Young Stars 



Our preliminary selection criterion of strong fractional NUV and FUV flux coinci des 
with the UV emission levels of the strong X-ray emitters of SLR09 and IRiaz et al.l (120061 ) as 
shown in Figure |3j This implies a rough upper limit to the age of our sample of 300 Myr as 
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estimated for the two ROSAT samples. After identifying the strongest UV emitters in our 
sample (F FUV /Fj > 1CT 5 and F NUV /Fj > 10~ 4 ), we used the same spectroscopic age-dating 
criteria discussed in detail in SLR09 and summarized in Figure El i.e. low gravity, lithium 
absorption and strong Ha emission as diagnostics of youth. 



6.1. Spectroscopic Youth Indicators 



Models of PMS stars show that lower-mass stars take longer to contract to the main- 
sequence (MS), e. g. a 0.5 Mq star wil l reach the MS within 100 Myr whereas a 0.1 M star will 
do so in 300 Myr ( IBaraffe et al.lll998[ ) and thus determining if a M star has l ow surface gravity 



provid es an upper limit to its age. We used the CaH gravity index from iKirkpatrick et al. 



(119911 ) defined as the ratio of the mean intensity in two passbands, a 'continuum' band and a 
molecular absorption band of CaH A6975: [7020-7050 A]/ [6960-6990 A]. Since we have 15-20 
times the resolution of previous M dwarf surveys, we also used a narrower 5-A CaH index 
[7044-7049]/ [69 72. 5-6977.5] providing a more discriminating scale with which to identify 
low-gravity stars (SLR09). Both indices are plotted as function of SpT in Figure [3 

An important caveat to using the TiO and CaH molecules as temperature and gravity 
diagnostics is their dependence on metallicity. Higher metallicity will mimic later spectral 
types and lower surface gravities. (See discussion in Secti on 5.1.1 of SLR09.) We thus 
also used the K I A7699 A line as a gravity indicator (e.g. ISlesnick et al.l 120061 ) . Care is 
required with this line as well as it is affected by stellar activity such that higher levels of 
chromospheric emission fill in the absorption cores and reduce the measured EWs. 

Combining the effects of the chromosphere on K I with the uncertainties in metallicity 
on the TiO and CaH indices, we considered a target as having low-g only if both the CaH 
and K I measurements indicate that it is so. We flagged a target as such in Table H] if it 
falls on or below (within error bars) the best-fit curves to the observed /3 Pic members in 
Figures [7] and [SI Data for the 12-Myr old Pic members were taken from SLR09. Out of 
the 16 UV-bright stars with Ha in emissio n, 8 have low surface gravity with upper age limits 
of less than 160 Myr (jBaraffe et al.l 119981 ). Three of these 8 show additional indications of 
youth and will be discussed in more detail in the following sections^ 



Lower limits on the stellar ages for early M dwarfs are provided for those stars with 



8 The 5 stars with no additional signs of youth have SpTs ranging from M0. 3 to M3.9. For these early 
SpTs, the difference in their predicted masses for a given age produce negligible results in the upper age 
limits set by the models. If the SpTs were known more precisely, the most significant difference would be in 
the upper age limit set for the M3.9 star, which would change from 160 Myr to 120 Myr. 
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no lithium absorptio n (A6708 A) using the lithium depletion time scales calculated by 
Chabrier et al.l (ll996[ ). However, it has been recently shown empirically for at least the 
P Pic moving group, that lower age limits of individual stars based on the lack of lithium 
absorption systemat ically over-estimates the star's age as comp a red t o model isochrones 
(jYee fc Jensen! |2010| ). In agreement with this. iBaraffe fc Chabrierl (120101 ) recently presented 
models where stars that are exposed to episodic accretion have higher internal temperatures 
and thus enhanced lithium depletion compared to stars of the same age and mass but with- 
out such accretion. This would imply that the stars with lower age limits in Table H] may 
indeed be younger, perhaps even as young as ~10 Myr. 

We have measured large lithium EWs^l in 2 of our targets (TWA candidates 1207-3230 
and 1226-3316; see Section [7]) and one marginal detection (0.18 ± 0.05 A) in candidate 
1131-4826 (Figure ED setting an upper age limit just from the lithium detection of ~40 Myr 
for the first two and 15 Myr for the last due to its earlier SpT. 

Lastly, the most stringent upper l imit is provided by the detection of active stellar accre- 



tion. IBarrado y Navascues fc Martini (120031 ) have produced an empirical accretion diagnostic 
by calibrating the Ha EW as a function of SpT. This accounts for the increase in EW simply 
due to the drop in continuum flux from cooler stars (Figure fTUl) . Though the accretion curve 
is not thought to be very robust for objects of SpT later than M5.5, due to the few late-M 
cluster members used to calibrate the sequence, it does serve as an outer envelope of the 
chromospheric emission in early Ms. In addition to the EW limits, IWhite fc Basril (120031 ) 
impose a Ha 10% velocity width criterion of >270 km s" 1 for optically veiled T Tauri stars, 
and >200 km s" 1 for non-optically veiled T Tauri stars. 

We plot the Ha EW of our target stars in Figure [10J One of the two strong lithium stars 
exhibits extremely large Ha emission, and is discussed in more detail below. Two additional 
stars in Figure [TU] have Ha EWs beyond this accretion/non-accretion boundary, both of 
which are SB2s, including TWA 3Aab of which at least one component is still accreting 
as determined by its very broad Ha velocity width (395 km s _1 ). The other SB2, TWA 
candidate 1013-3542, must have enhanced chromospheric emission due to the tidal spin-up 
of the two stars, not accretion, as there are no additional spectral signatures of youth, i.e. no 
lithium absorption nor any indication of low gravity. 



9 The lithium abundances have not been corrected for possible contamination with the Fe I line at 6707.44 
A. Uncertainties in the setting of continuum levels prior to measurement induce EW errors of about 0.01-0.02 
A with a dependence on the S/N in the region. We therefore consider our 2a detection limit to be 0.1 A. 
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6.2. Kinematics of the TWA Candidates 



The high resolution of the data provided RV measurements to better than 1 km s _1 in al- 
most all cases, which we used in conjunction with the star's photometr ic distance and proper 
motio ns (Figure [TTj) to measure its 3-dimensional space velocity (UVW; I Johnson fc Soderblom 
19871 ). This provides a promising way to determine stellar ages by linking stars kinematically 
to one of the several known YMGs or associations, including TWA. 

We calculated photometric distances for our sample of TWA candidates, as well as for 
the known TWA members we observed and have listed them in Tables [3] and [51 We used the 
2MASS K magnitud e, R — K or V — K ( when available) colors and the age and absolute K 



magnitudes from the iBaraffe et al.l ( 119981 ) models to measure the photometric distance. We 
took into account the spectroscopically determined ages (or upper age limits) from Table IH 
since PMS stars are over-luminous compared to dwarfs and will appear closer than they are. 
For the proposed two new TWA members (discussed belo w), we adopted the 10 Myr age 
of TWA (IBarrado Y Navascuesl 120061 ; iMentuch et al.l 120081 ). Distances are corrected for any 
known unresolved binaries. Given the uncertainties in the models, the metallicities of the 
stars, and the age ranges provided in Table El we estimate the errors to be ^20% for the 
non-TWA members and ~10% for the TWA members, which have more precise ages, and 
not only upper age limits. 



The UVW velocities for the targets are shown in Figure [121 Seven TWA candidates 
fall near the 2a error ellipse for the UVWs of known TWA members. Their RA/DEC co- 
ordinates are 1037-3505, 1039-3534 A, 1039-3534 B, 1130-4628, 1131-4826, 1207-3230, and 
1226-3316. And given the large uncertainty in distance, we plot the UVWs of these 7 with 
a range of possible distances (30-130pc) in Figure [121 

All 7 of the stars appear to have low surface gravity, with only the last two exhibiting 
strong lithium absorption. (See Section [3) The eighth star listed as low-g in Table HI 
2MASS1111-3937, has strong and broad Ha emission as w ell but UVW velocities inco nsistent 
with the "good UVW box" defined for young stars by IZuckerman fc Songl (120041 ) . It is 
thus probable that 2MASS1111-3937 is an unresolved SB2 with broadened spectral features 
caught at an orbital phase near conjunction. 

Candidate 1131-4826 has a weak lithium detection with an EW of 0.18 ± 0.05 A, setting 
an age limit of 15 Myr. Based on its age, distance of 93 ± 19 pc, RV of 17.02 ±1.16 km s -1 
and UVW velocities (-9.0 ± 3.3, -22.6 ± 19, -4 1 ± 1.8 km s" 1 ), we conclude that this 
target is a member of LCC ( jde Zeeuw et al.lll999l ; iMamajek et al.l 120021 : iBitner et al.l 120 101 ) 
rather than TWA. It is also worth pointing out that LCC 1131-4826 appears to still be 
accreting based on its broad Ha profile (233 km s _1 ; Figure IT4"|) . making it the first known 



accreting M star in LCC (IPreibisch fc Mamajekll2008l ). 



The remaining four low-g stars with UVWs near the TWA UVW error ellipse show no 
additional signs of youth beyond the UV excess and Ha emission, and have upper age limits 
ranging from 110 - 300 Myr. They range in distance of 41 - 84 pc and, assuming lithium is 
a necessary youth indicator, are not obviously part of either TWA or LCC. 



Two New TWA Members: TWA 31 and TWA 32 



Two of our targets share all the same spectroscopic, photometric and kinematic charac- 
teristics of known TWA members including low surface gravity, strong Li absorption, strong 
Ho; emission, plus RVs and UVWs consistent with previously known TWA members. These 
two likely members have 2 MASS coordinates [12:07:10.89 -32:30:53.72] and [12:26:51.35 - 
33:16:12.47], SpTs of M4.2 and M6.3, and in keeping with tradition, we dubbed them TWA 
31 and 32, respectively. They are identified in Tables [T] and HJ and are marke d by large circles 



in th e figures. Proper motions for the two are from the NOMAD catalog (IZacharias et al. 



20051 ). and agree well with the proper motion vectors of known members (Figure [TT 



The average Li EW of undisputed members TWA 1-12 from (IMentuch et al.ll2008l ) is 
0.52 with rms=0.06 A, and TWA 32's EW is consistent with this (0.60 ± 0.05 A). TWA 
31 has a sligh tly lower-than-average Li EW (0.41 ± 0.05 A), likely due to optical veiling 
(jDuncanlll99ll ). TWA 31 also has by far the strongest Ha emission in our sample with an 
Ha EW of -115 A and an extremely accretion-broadened 10%- velocity width of 447 km s _1 
(Figure [To^ . characteristics comparable to TW Hydrae itself We conclude that TWA 31 
is also an accreting T Tauri star with an age of <10 Myr. TWA 31 also emits strongly 
at He I ( EW of -3.6 A at A66 78 A and -10.3 A at A5867 A), yet another indication of 
accretion ( IMohanty et al.ll2005l ). making it only the 6th known TWA a ccretor - the other 
ones being, TW Hyd Hen 3-600, TWA 14 flMuzerolle et~al1l2000bl . l20~0~lT ) and TWA 30A+B 
f lLooper et aUl2010al lrl. 



The photometric distances for TWA 31 and 32 are 110 ± 11 pc and 53 ± 5 pc (taking 
binarity into account; see below), respectively. Although the distance to TWA 31 is relatively 
large compared to most of the known TWA members, it does not appear to be part of 
any other neighboring association and is certainly too young to be a member of the 120- 



'Value s for TW Hyd are: \og(F NUV / Fj) = - 2.081 log( F FUV / Fj) = -2.456, Li EW = 0.467 ± 
0.021 A jMentuch et al1 l2008l). Ha EW = 220 A jlleidl 120031) . Ha 10% velocity width = 400 km s" 1 
(jAlencar fc Batalhalbopj ). 
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pc, 16-Myr old LCC, which is adjacent on the sky to TWA. It is possible that the large 
distance of TWA 31 implies that it is part of an unidentified young star association in 
the direct ion of TWA, as specula ted about several other distant, yet clearly young, TWA 



members ( [Weinberger et al 



20111 ) . 



7.1. Keck LGS AO Imaging of TWA 32 



We imaged TWA 32 on UT 2010 May 22 using the sodium laser guide star adap- 
tive optics (LGS AO) s y stem of the 10-meter Keck II Telescope on Mauna Kea, Hawaii 
( jWizinowich et al.l 120061 ; Ivan Dam et al.l 120061 ) . We used the facility IR camera NIRC2 
with its 10.2" x 10.2" field of view during photometric conditions. The LGS provided the 
wavefront reference source for AO correction, with the tip-tilt motion measured simulta- 
neously from the star itself. We obtained a set of dithered images wi th the broadband K 
(2.20 p) filter from the M auna Kea Observatories filter consortium ( iSimons fc Tokunaga 
20021 ; iTokunaga et al.l 120021 ) and easily resolved the target into a nearly equal-flux binary. 
Images were reduced in a standard fashion, and the relative ast rometry and pho tometry were 



derived using a multi-gaussian representation of the PSF (e.g. iLiu et al.l 120081 ) . Astrometry 
was corrected fo r instrumen t al dis tortion, with the absolute calibration of plate scale and 
orientation from iGhez et al.l ( 120081 ). We measured a separation of 656.1 ± 0.4 mas, a posi- 
tion angle of 11.65 ± 0.08 degs, and a flux ratio of 0.230 ± 0.008 mag, with the uncertainties 
determined from the scatter in the individual images and the overall astrometric calibration 
uncertainties. 



8. Summary 

We set out to find the 'missing' mid-M dwarfs in TWA by cross-matching optical (HST 
GSC), infrared (2MASS) and UV (GALEX) catalogs, filling in the stellar mass function 
in the association and providing excellent targets for direct imaging searches for substellar 
companions and circumstellar disks. GALEX provides a new and more sensitive resource that 
enables a major expansion of the young low-mass census, far beyond previous data sets. We 
found that NUV data yield many candidates, but applying FUV criteria provide an excellent 
means to distinguish between the (never-before-delineated) quiescent emission of old stars 
{Ffuv/Fj < 10 -5 ), the faint sources (FUV not detected), and the truly young targets with 
high levels of UV emission (F FUV /Fj > 10~ 5 and F NUV /Fj > 10" 4 for R - J > 1.5). The 
photometric cross-matching yielded 34 UV-bright low-mass stars with SpT between M0 and 
M6 within ^100 pc and 1000 sq. deg. of the TWA with ages likely less than 300 Myr. 



-14- 



Ground-based optical low-resolution spectroscopy of 30 identified 16 with Ha emission 
which were followed up with high-resolution spectroscopy!"! Of these, 2 (possibly 3) are 
old SB2s with tidally enhanced UV emission. Six are nearby field Ms with ages probably 
younger than 300 Myr, based on their strong UV emission with no additional signs of youth. 
One candidate appears to be an accreting new member of the 16-Myr old LCC, and 5 are 
low-gravity M dwarfs with maximum ages ranging from 110 - 300 Myr. Four of these 5 
low-g stars are kinematically identical to the previously known TWA members, yet are likely 
older than 20 Myr based on the absence of Li absorption. Thus identifying new YMG 
members based on kinematics and strong X-ray or UV emission alone may not be sufficient 
and spec troscopic observations are n ecessary for confirmation. However, in light of the recent 



work by iBaraffe fc Chabrierl (120101 ). these stars may be much younger and possibly TWA 



members despite not having any Li absorption. 

Lastly, 2 stars in our sample exhibit all the spectroscopic, photometric and kinematic 
characteristics of ~10 Myr-old TWA members including low surface gravity, strong Li absorp- 
tion, strong Ha emission and UVW velocities. These new members, TWA 31 (SpT=M4.2) 
and TWA 32 (SpT=M6.3), have photometric distance of 110 ± 11 pc and 53 ± 5 pc, re- 
spectively. Followup Keck/LGS AO observations resolved TWA 32 into two near equal-flux 
(0.230 ± 0.008 mag in K) stars with a separation of 656.1 ± 0.4 mas. TWA 31 also exhibits 
an extremely accretion-broadened Ha profile (447 ± 10 km s ~ l ) with a slig htly lower-than- 



average Li EW (0.41 ± 0.05 A), likely due to optical veiling ( )Duncanlll99lf ). making it only 
the 6th known active accretor in TWA. 

Our new GALEX/ AIS search method successfully recovered 2/3 of the known TWA 
members (corresponding to the 2/3 of the sky covered in the GR4 data archive release), 
making it surprising that only 2 new mid-M members were discovered. With the peak of the 



spectral-type (and mass) function at SpT=M3-M4 (IBochanski et all 120101 and references 
therein), the expected number of newly found M dwarfs is substantially greater than 3, 
likely closer to 20 (based on the known number of early-Ms in TWA). Our results imply that 
either TWA has an unexpected mass function, or a significant fraction of 10-Myr M dwarfs 
have depleted all their lithium and were eliminated from the membership list. This latter 
possibility would imply that low-g stars that are kinematically identical to TWA but lacking 
Li may indeed be bona fide members of the association. 

E.S. thanks Bernie Shaio and Tony Rogers of MAST/STScI for GALEX/ AIS query 
support. Also, we appreciate the helpful comments on the manuscript by the referee, Eric 



L Several of the UV bright targets with no Ha emission are white dwarf+M dwarf pairs. 
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All Sky Survey, with access to the last two provided by Vizier and SIMBAD. 2MASS is a 
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u.oyo 
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ZO.OOo 
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q nn 
-o.UU 


no 
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OSO 8QA 
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U.ZOl 
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0S7 AC\A 
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07 SAO 
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1Z.00 
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n i oo 
u. iyo 


U.UOo 
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- 1 .oo 


Q fiO 
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no 
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OSO /inA 

zoy .4uo 


OQ OQA 

zo.zoo 


1 S K7 
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1 A Ol Q 
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U.0Z1 


U.O40 


OIQ 
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1 OQ 
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10 nn 07 1^1 Q/i nt; Q7 17 

1Z UU Z/.01 -o4 UO oi.li 


OOI flQC, 

zyi.uoo 


07 t;ns 

z / .oyo 


i q Qf; 
lo.oO 


S 70Q 

o. / Zo 


n koo 

U.OZZ 


n oai 

U.ZOl 


Q SA 

-o.oO 


A RA 
-4.04 


no 




12 03 08.07 -38 26 55.54 


292.643 


23.465 


11.66 


8.565 


0.594 


0.271 


-3.77 


-4.54 


yes 


ROSAT source 32" away 


12 07 10.89 -32 30 53.72 


292.208 


29.457 


16.56 


12.115 


0.558 


0.375 


-2.89 


-2.81 


yes 


TWA 31 


12 19 07.68 -41 01 57.81 


296.393 


21.441 


16.93 


13.012 


0.566 


0.303 


-1.74 


-2.42 


yes 




12 19 29.82 -34 14 24.54 


295.447 


28.178 


19.07 


14.347 


0.711 


0.331 


-1.69 


-2.15 


no 




12 20 25.78 -43 04 06.97 


296.934 


19.442 


17.97 


13.719 


0.589 


0.367 


-1.99 


-2.12 


no 




12 21 33.27 -41 40 29.12 


296.968 


20.865 


17.97 


13.866 


0.706 


0.261 


-1.43 


-1.35 


no 




12 24 08.36 -31 34 27.48 


296.156 


30.940 


18.11 


13.850 


0.624 


0.250 


-2.57 


-2.50 


no 




12 26 43.75 -41 47 37.31 


298.004 


20.856 


12.54 


9.082 


0.817 


0.267 


-1.55 


-3.63 


no 




12 26 51.35 -33 16 12.47 


297.040 


29.320 


14.92 


9.783 


0.569 


0.339 


-3.38 


-3.32 


yes 


TWA 32 


12 33 43.57 -32 51 26.29 


298.636 


29.877 


16.95 


12.149 


0.556 


0.317 


-2.62 


-2.41 


yes 




12 45 43.60 -37 04 36.24 


301.673 


25.790 


18.85 


13.649 


0.556 


0.398 


-2.14 


-2.40 


yes 




12 56 49.64 -30 07 37.37 


304.317 


32.730 


16.95 


12.883 


0.655 


0.294 


-2.60 


-2.82 


no 





a Note that Fj is calculated using the full width of the 2MASS J band filter, 0.3 /im. 



Table 2. TWA members observed with MIKE 



Name 


TP I TP \ 

logyr NUV 1 r J ) 


TP 1 TP \ 

logyi 1 FUV / r J) 


op 1 

l\ I 

M 


CaH-narr. 
Index 


T T?"\Xf 

X 

A 


T ; T?\X7 

X 

A 


TJrt, T?"\Xf 

X 

A 


Binarity 


TWA 2AB 


-3.87 


-4.62 


1.5 


1.21 


0.69 


0.52 


-1.72 


VB (Brandeker et al. 2003) 


TWA 3Aab 


-3.194 a 


-3.64 


3.9 


1.31 


0.87 


0.51 


-40.89 


VB. SB2 (Muzerolle et al. 2000a) 


TWA 3B 


-3.194 a 


-3.64 


3.9 


1.30 


0.85 


0.54 


-4.26 


VB 


TWA 5Aab 


-3.57 


-4.19 


1.9 


1.33 


0.90 


0.63 


-6.37 


VB. SB2 (Torres et al. 20061 


TWA 7 


not observed 




2.4 


1.28 


0.81 


0.55 


-5.39 




TWA 8A 


-3.61 


-4.36 


2.4 


1.37 


0.89 


0.55 


-5.04 


VB 


TWA 8B 


-3.72 




5 b 


1.29 


0.85 


0.58 


-6.21 


VB 


TWA 10 


-3.87 


-4.50 


2.6 


1.36 


0.89 


0.50 


-5.46 




TWA 12 


-3.72 


-4.39 


1.6 


1.21 


0.83 


0.53 


-5.10 




TWA 14ab 


-2.90 


-3.42 


0.6 


1.21 


0.84 


0.59 


-5.68 


SB2 (Jayawardhana et al. 2006) 


TWA 15B 


not observed 




2.2 


1.38 


1.06 


0.55 


-9.64 


VB 


TWA 16 


-3.81 


-4.18 


1.8 


1.32 


0.78 


0.38 


-3.08 




TWA 22AB 


not observed 




6.5 


1.52 


1.69 


0.65 


-10.48 


VB fBonnefov et al. 2009) 


TWA 23ab 


-4.05 


-4.65 


2.9 C 


1.29 


0.70 


0.50 


8.12 


SB2 (this work) 



a G ALEX cannot resolve TWA 3 A and 3B. 

b The TiO index gave SpT=M3, but we list the published SpT from lTorres et al.1 j2003h . 
c Both components of TWA 23 have SpT=M3. 



Table 3. Kinematics of known TWA members observed with MIKE 



Name 


RA & DEC 


pmRA a 


pmDec 


Dist. b 


RV 


U 


V 






W 










dej 




mas yr _1 


mas yr 


pc 


km s~ 1 


km s~ 1 


km 


s 






km 


s 






TWA 


2AB 


167.31 


-30.03 


-yo.o it z.y 


-ZO.0 It Z.o 


To 


lU.oo it U.ol 


1 A Q -I- 1 7 
-14. y ± l.i 


-18.2 


± 


1 

1 


1 

. 1 


-7.7 


± 


1 

1 . 


1 

. 1 


1 WA 


3Aab 


167.62 


-37.53 


i no -t- 7 

-1UU It 1 


14-1-11 
-14 It 11 


al 


Q M 1 (1 Of! 

y.oz it u.oo 


10 4-1 7 
- 1U.U It 1. / 


-14.0 


± 


1 
1 


1 

. 1 


-3.9 


± 


1 

1 


.0 


TWA 


3B 


167.62 


-37.53 


-100 ± 7 


-14 ± 11 


31 


9.89 ± 0.62 


-9.9 ± 1.6 


-14.4 


± 





.9 


-3.8 


± 


1 


.6 


TWA 


5Aab 


172.98 


-34.61 


-85.3 ± 3.6 


-23.3 ± 3.7 


38 


13.30 ± 2.00 


-8.5 ± 1.4 


-18.7 


± 


1 


.9 


-2.5 


± 


1 


.3 


TWA 


7 


160.63 


-33.67 


-122.2 ± 2.2 


-29.3 ± 2.2 


34 


12.21 ± 0.24 


-13.2 ± 1.5 


-17.7 


± 





.7 


-8.5 


± 


1 


.4 


TWA 


8A 


173.17 


-26.87 


-90 ± 2 


-20 ± 15 


14 


8.34 ± 0.48 


-13.1 ± 2.1 


-15.8 


± 


1 


.6 


-4.4 


± 


2 


.6 


TWA 


SB 


173.17 


-26.87 


-86 ± 3 


-22 ± 38 


27 


8.93 ± 0.27 


-6.8 ± 2.4 


-12.7 


± 


2 


.1 


-0.7 


± 


3 


.9 


TWA 


10 


188.77 


-41.61 


-78 ± 23 


-32 ± 8 


67 


6.75 ± 0.40 


-15.7 ± 6.6 


-21.1 


± 


4. 


.2 


-8.6 


± 


2 


.7 


TWA 


12 


170.27 


-38.75 


-60 ± 14 


-12 ± 25 


63 


13.12 ± 1.59 


-11.0 ± 5.1 


-19.1 


± 


2 


.8 


-4.6 


± 


6 


.7 


TWA 


14ab 


168.36 


-45.40 


-43.3 ± 2.6 


-7 ± 2.4 


113 


15.83 ± 2.00 


-14.4 ± 2.3 


-23.0 


± 


2 


.1 


-8.1 


± 


1 


.8 


TWA 


15B 


188.59 


-48.26 






100 


10.03 ± 1.66 




















TWA 


16 


188.73 


-45.64 


-53.2 ± 5.2 


-19 ± 5.2 


65 


9.01 ± 0.42 


-8.5 ± 1.9 


-17.2 


± 


1 


.4 


-4.0 


± 


1 


.7 


TWA 


22AB 


154.36 


-53.91 


-174.8 ± 9 


-13.6 ± 9 


20 


13.57 ± 0.26 


-10.1 ± 1.5 


-16.3 


± 





.1 


-9.7 


± 


1 


.3 


TWA 


23ab 


181.86 


-32.78 


-44 ± 7 


-12 ± 3 


61 


8.52 ± 1.20 


-7.0 ± 2.1 


-14.0 


± 


1 


.6 


-1.0 


± 


1 


.1 



a Proper motions are from the NOMAD catalog ( Zacharias et al]|2005l) . 

b Ph otometric distances are calculated using the iBaraffe et al . (1998) models with an age of 10 Myr, T e ff from lMentuch et al] 
(2008) and taking into account bina rity assuming equ al flux components. Uncertainties are w 10%. Distances agree within error 
bars with the predicted values from Mamaiek (2005) and E. E. Mamajek (private communication), except TWA 15B which is 
predicted to be at 41 ± 6 pc, respectively. 



Table 4. Spectral Age Diagnostics of TWA candidates observed with MIKE 



RA & 


DEC 




SpT 


CaH-narrow 


K I EW 


Li EW 


Ha EW 


Ha 10% width 


Youth 


Age b 


Note 


2MASS 






M- 


Index 


A 


A 


A 


km s _1 


Indcx a 


Myr 












(±0.5) 


(±0.03) 


(±0.2) 


(±0.05) 


(±0.5) 


(±4) 








10 13 21.43 


-35 


42 


oo.yo 




1 


.38 


z.oo 


<U.l 


OK Q 
-ZO.O 


Z ((J 




old 


SB2, rapidly orbiting 


10 37 10.47 


-35 


05 


ni c;o 
Ul.OZ 


1 

o.Z 


1 


.40 


1 /lO 

i .^y 


<U.l 


-0. 1 


lUo 


1 nnn 
1UUU 


iO — oUU 




10 39 52.76 


-35 


34 


Uo.UvS 


n f : ; 
U.u 


1 


.18 


n eo 

U.oZ 


^ ft 1 


1 

-1.0 


o 1 


1 1 nn 


on i in 
zu — 11U 


V Jj (IN ) 


10 39 52.76 


-35 


34 


HQ HQ 

Uo.Uo 


u . o 


1 


.16 


U.oo 


<U.l 


1 % 
-1.0 




1 1 nn 
11UU 


on i in 
zu — 11U 


V rS [o ) 


11 02 53.73 


-31 


45 


10.57 


2.3 


1 


.38 


1.67 


<0.1 


-3.2 


79 


0000 


20 - 300 




11 11 46.37 


-39 


37 


34.78 


2.2 


1 


.52 


1.09 


<0.1 


-12.4 


124 


1101 


25 - 130 


SB2? C 


11 11 52.67 


-44 


01 


53.87 


3.9 


1 


.42 


1.13 


<0.1 


-4.8 


91 


1100 


90 - 160 




11 18 12.37 


-32 


35 


59.09 


4.2 


1 


.46 


2.09 


<0.1 


-5.0 


79 


1000 


90 - 300 




11 30 53.55 


-46 


28 


25.19 


2.1 


1 


.37 


1.31 


<0.1 


-1.1 


77 


1100 


20 - 130 




11 31 14.83 


-18 


26 


27.98 


3.5 


1 


.35 


1.18 


0.18 


-7.3 


233 


1111 


«15 


LCC member 


12 03 08.07 


-38 


26 


55.54 


0.7 


1 


.30 


0.61, 0.26 


<0.1 


-1.9 


95, 73 




>300 


SB2 


12 07 10.89 


-32 


30 


53.72 


4.2 


1 


.35 


1.35 


0.41 d 


-114.8 


447 


1111 


10 


TWA 31 


12 19 07.68 


-11 


01 


57.81 


2.9 


1 


.35 


1.46 


<0.1 


-3.3 


72 


1000 


40 - 300 




12 26 51.35 


-33 


16 


12.47 


6.3 


1 


.39 


1.21 


0.60 


-12.6 


127 


1111 


10 


VB, TWA 32 


12 33 43.57 


-32 


51 


26.29 


4.6 


1 


.51 


2.54 


<0.1 


-2.6 


73 


0000 


180 - 300 




12 45 43.60 


-37 


04 


36.24 


5.5 


1 


.66 


3.07 


<0.1 


-4.1 


144 


0000 


180 - 300 





a Youth index in binary format in order of least to most restrictive age indicator: low-g from CaH, low-g from K i, Li detection, accretion-level Ha emission. 

b Lower li mits on the stellar ages for early M dwarfs are provided for those stars with no lithium absorption (A6708 A) using the lithium depletion time scales 
calculated bv lChabrier et ah | lll996h . However, it has been recently shown empirically, for at least the 12-Myr old /3 Pic moving group, tha t lower ages lim its of 
individual stars based on the the lack of lithium absorption systematically over-estimates the star's age as compared to model isochrones jYee &: Jensenlfeoiol) . 
This would imply that the lower age limits may be even lower, perhaps even as low as the accretion limit of fslO Myr. The upper age limit is set by the UV 
emission and/or low gravity using the evolution models of lBaraffe et al.l jl998h . 

c This star may be an unresolved SB2. See Section 16.21 for more details. 

The lower value of Li EW compared to other TWA members is likely due to veiling due to accretion touncanlll99lT) . 



Table 5. Kinematics of TWA candidates observed with MIKE 







RA & 


DEC 






pmfxA 


pmDec 


"DTP TP a 


uisi. 


XXV 


I T 

u 


V 


VV 


Note 






2MASS 








mas yr 


mas yr 




pc 


km S 


km s 


km S 


km s ^ 




10 


13 


21.43 


-35 


42 


36 


.93 


18 ± 6 


-14 ± 11 


1 


297 










SB2 


10 


37 


10.47 


-35 


05 


01 


.52 


-52.8 ± 2.4 


8.9 ± 2.4 


2 


41 


12.45 ± 0.54 


-8.7 ± 1.9 


-13.7 ± 0.7 


0.6 ± 0.7 




10 


39 


52.76 


-35 


34 


03 


.03 


-58.3 ± 1.3 


5.1 ± 4.2 


2 


82 


12.45 ± 0.47 


-19.1 ± 4.0 


-16.8 ± 1.2 


-5.1 ± 2.3 


VB (N) 


10 


39 


52.76 


-35 




03 


.03 


-58.3 ± 1.3 


5.1 ± 4.2 


2 


82 


12.82 ± 0.21 


-19.1 ± 4.0 


-17.1 ± 1.1 


-5.0 ± 2.3 


VB (S) 


11 


02 


53.73 


-31 


45 


10 


.57 


22 ± 3 


4 ± 4 


1 


218 


95.27 ± 0.22 


28.6 ± 4.9 


-76.8 ± 2.4 


53.9 ± 4.4 




11 


11 


46.37 


-39 


37 


34 


.78 


-34 ± 12 


2 ± 12 


1 


200 


-72.00 ± 1.24 


-43.8 ± 11.2 


55.9 ± 4.6 


-34.3 ± 10.7 


SB2? 


11 


11 


52.67 


-11 


01 


53 


.87 


-22 ± 2 


-12 ± 4 


1 


34 


17.64 ± 0.30 


2.1 ± 0.6 


-17.9 ± 0.4 


1.6 ± 0.9 




11 


18 


12.37 


-32 


35 


59 


.09 


-38 ± 3 


8 ± 8 


1 


106 


-0.70 ± 0.66 


-18.6 ± 4.3 


-4.7 ± 1.9 


-3.7 ± 3.4 




11 


30 


53.55 


-46 


28 


25 


.19 


-33.7 ± 3.2 


1.1 ± 1.8 


2 


62 


10.03 ± 0.12 


-5.6 ± 1.9 


-12.9 ± 0.8 


-0.3 ± 0.8 




11 


31 


14.83 


-48 


26 


27 


.98 


-40.2 ± 3.3 


-5.8 ± 1.5 


2 


93 


17.02 ± 1.16 


-9.0 ± 3.3 


-22.6 ± 1.9 


-4.1 ± 1.8 


LCC 


12 


03 


08.07 


-38 


26 


55 


.54 


106.6 ± 1.8 


-18.3 ± 1.1 


2 


45 


-49.64 ± 1.00 


3.6 ± 4.3 


51.2 ± 2.0 


-19.0 ± 0.5 


SB2 


12 


07 


10.89 


-32 


30 


53 


.72 


-42 ± 6 


-36 ± 3 


1 


110 


10.47 ± 0.41 


-8.8 ± 3.2 


-25.5 ± 2.5 


-14.6 ± 3.4 


TWA 31 


12 


19 


07.68 


-41 


01 


57 


.81 


-74 ± 20 


34 ± 4 


1 


101 


6.44 ± 0.11 


-32.6 ± 10.9 


-18.2 ± 5.3 


12.9 ± 3.0 




12 


26 


51.35 


-33 


16 


12 


.17 


-62.2 ± 3.5 


-24.7 ± 3.9 


2 


53 


7.15 ± 0.26 


-8.6 ± 1.4 


-15.7 ± 1.1 


-3.4 ± 1.1 


VB, TWA 32 


12 


33 


43.57 


-32 


51 


26 


.29 


-18.4 ± 5.8 


31.5 ± 5.8 


3 


34 


10.38 ± 0.29 


0.2 ± 1.2 


-7.4 ± 0.7 


9.4 ± 1.2 




12 


45 


43.60 


-37 


04 


36 


.24 


21.7 ± 9.6 


-80.7 ± 9.8 


3 


73 


-10.81 ± 1.82 


8.2 ± 4.1 


2.3 ± 2.8 


-29.7 ± 6.0 





a Proper motion references: 1 = IZacharias et al~ll2005l . 2 = IZacharias et al]|2010l . 3 = lRoeser et aLlboid 

b Photometric distances take into account youth (using Baraffc e t al.|[l998h and binarity assuming equal flux components. Uncertainties are ss 20% for 
the non-TWA members and «10% for the TWA and LCC members. See text for more details. 
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Fig. 1. — Spectral type histogram of known TWA members including TWA 31 and 32. 
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Fig. 2. — A 2MASS color-color plot of our TWA candidates plus known TWA members 
The two stars circled are TWA 31 and TWA 32, and are discussed in detai l in Section CTJ 
Neithe r shows significant K-band excess using the M dwarf color sequence of ICushing et al. 
J2005h . 
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Fig. 3. — Fractional near-UV flux plotted a gainst R — J for t he GALEX-detected M dwarfs 



including those from the NStars 25-pc list ( Reid et al.l 120071). The red squares represent X- 
ray-bright stars within 50 pc (< 300 Myr; Ishkolnik et al.lbood : iRiaz et al.ll20n3 ). The TWA 
candidates targeted for spectroscopic followup have both high NUV and FUV fractional 
luminosities. However those with Fnijv/ Fj > 0.01 with no H a emission are most likely M 
dwarf + white dwarf pairs (green crosses; ISilvestri et al.ll2007l ). The two new TWA members 
reported here are identified with large circles. 
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Fig. 4. — Galactic latitude and longitude of known TWA members and our candidates. 
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Fig. 5. — A color-magnitude diagram showing the TWA candidates for which we collected 
low- resolution data. The lines are the applied color and magnitude cuts detailed in Section [3] 
which set a photometric distance limit of 100 pc (assuming a sample of main-sequence stars). 
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Fig. 6. — A summary of age diagnostics used in this paper. Each technique provides an 
upper limit and in the case of lithium, a lo wer limit if none is detected. The limits set by 
low-gravit y are from evolutionar y models oflBaraffe et al.l (1998 ) a nd lith ium depletion from 
models of Chabrier et al. (jl996 ). Barrado y Navascues fc Martin ( 2003 ) set an upper limit 
of 10 Myr for a star still undergoing accretion. 
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Fig. 7. — The CaH indices (wide on the left; iKirkpatrick et al.l (jl99ll ). and narrow on the 
right; SLR09) for our sample of stars with high resolution optical spectroscopy. The red 
dashed curves are polynomial fits to our CFHT and Keck observations of (3 Pic M dwarfs 
(SLR09): CaR wide = -0.0067 SpT 2 + 0.0986 SpT + 1.0633 and CaH ^ = -0.00891 SpT 2 
+ 0.13364 SpT + 1.1053. The rj Cha data are from lLvo et aJ fho04f ). 
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Fig. 8. — K I equivalent widths as a function of spectral type. The red dashed curve is a 
polynomial fit to the (3 Pic observations from CFHT and Keck (SLR09): EW K i = -0.02821 
SpT 2 + 0.20731 SpT + 0.62911. 
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Fig. 9.— Lithium EWs of TWA 31, TWA 3 2 and LCC 1131-4826 compared with those of 
known TWA members (IMentuch et al.ll2008l ). 



-120 
-100 

_ -80 

»< 

£ -60 
W 

.3 -40 
-20 




i 1 i 1 r 

■ |3Pic 

O TWA candidates 
* TWA members 



SB2 O 



SB2 * 



_l i I i I i I i I i L 

1 2 3 4 5 
SpT (M subclass) 



Fig. 10. — Ha equivalent widths as a function o f spectral type. The dashed curve repr esents 
the empirical accretion boundary determined by lBarrado y Navascues fc Martini (120031 ). The 
2 SB2s above the accretion curve are TWA 3Aab, which has at least one component still 
accreting, while the SB2 from our candidate list is not. 
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Fig. 11. — RA and DEC positions of candidates (left) and known TWA members (right) 
with proper motion vectors. Proper motions with references are listed in Table 16.21 
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Fig. 12.— Left: UVW velocities of the TWA candidates observed with MIKE. Right: A 
zoomed-in view of the 7 candidates clustered around TWA's UVW error ellipse (±2<r) cen- 
tered on the average UVW=-10.5, -16.9, -4.8 km s _1 . UVWs of the known TWA members 
listed in Table [3] are shown as asterisks. 
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Fig. 13. — UVWs calculated for a range of stellar distances for the 6 TWA candidates that 
fall near TWA's 2-cr error ellipse at their photometric distances. The distances start at 30 
pc (top right) and increase by 5 pc increments to a maximum of 130 pc. 
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Fig. 15. — Ha profiles of TWA 31. Three consecutive 900-s exposures were taken of TWA 
31 on UT 2009 June 07. The first observation in the series is the strongest emission profile 
followed by thick black and thin dashed curves. 



